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Abstract. The purpose of these notes is to provide the details of the Jacobian 
ring computations carried out in [l], based on the computer algebra system 
Magma g]. 

1. Finding an appropriate matrix 

Our first target is to construct the Jacobian ring associated to a Calabi-Yau X as 
described in section 2 of pQ. This variety X is defined by a matrix A £ C 4x8 subject 
to the condition that all 4 x 4-minors are non-zero. In practice, we work over the 
field Q or F p , p prime, instead of C. The matrices contained in the hyperelliptic 
locus were of special importance for us: These are of the form A := (a^) with 
fly = A*-, where (Aj., As) denotes a tuple of 8 distinct numbers in C. 

The Magma program below produces such an admissible matrix. There are three 
options: 

(i) Work with a user-defined matrix. In this case one sets randmat : =f alse 
and hyperell : =f alse. The line A:=RMatrixSpace(K,4,8) !0; (marked 
with "modify here") has to be replaced by a statement that defines the 
desired matrix, e.g. 

A:=Matrix(K, [ 

[ 1, 1, 1, 1, 1, 1, 1, 1], 

[ 1, 2, 3, 4, 5, 6, 7, 8], 

[ 1, 4, 9, 16, 25, 36, 49, 64], 

[ 1, 8, 27, 64,125,216,343,512]]); 



In this case the non-degeneracy condition is not checked; it has to be 
ensured by the user. 

(ii) Use a random matrix. Set randmat : =true and hyperell : =f alse. The pro- 
gram creates a random matrix in K 4x8 which satisfies the non-degeneracy 
condition. The parameter randrange specifies the maximal absolute value 
of the matrix coefficients that are created. 

(iii) Work with a user-defined matrix contained in the hyperelliptic locus. In 
this case put hyperell : =true, the value of randmat is irrelevant. In the 
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line marked with "modify here" the definition la:= can be replaced by any 
tuple of 8 distinct elements in K. 

For our implementation we used MAGMA version 2.11-2. It was carried out under 
the operating system Linux on an ordinary personal computer. The following 
program produces the desired matrix. 

randmat := true; 

randrange : = 10 ; 
hyperell : = true ; 

K := RationalsO; 

A := RMatrixSpace(K,4,8) !0; // <== modify here 

genpos:=not (randmat) and not (hyperell) ; 
while not(genpos) do 

N : =randrange ; 

if hyperell then 
la: = [] ; 

for i:=l to 8 do 

la[i] :=Random(2*N-2)-N+l; 
end for; 

la:=[l,2,3,4,5,6,7,8] ; // <== modify here 
A:-[]j 

for i:=0 to 3 do 

A[i+1] : = [] ; 

for j : =1 to 8 do 

A[i+1] [j] :=la[j]-i; 

end for; 
end for; 

else 

A: = [] ; 

for i:=l to 4 do 
A[i] : = [] ; 
for j : =1 to 8 do 

A[i] [j] :=Random(2*N-2)-N+l; 
end for; 
end for; 
end if ; 

A:=Matrix(K,A) ; 
seq: = [l,2,3,4,5,6,7,8] ; 
set : =SequenceToSet (seq) ; 
subsets :=Subsets (set, 4) ; 
res : =true ; 

for su in subsets do 

seq: =SetToSequence (su) ; 
B:-[]j 

for i:=l to 4 do 
B[i] : = [] ; 
for j : =1 to 4 do 

B[i] [j] :=A[j] [seq[i]] ; 
end for; 
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end for; 
B:=Matrix(K,B) ; 

res:=res and not (Determinant (B) eq 0); 
end for; 
genpos :=res ; 
end while; 

In the first part of the program a random or hyperelliptic matrix is created according 
to the user specification. The second part is to check the non-degeneracy condition. 
This is done in the straightforward way: A loop runs through all 4-element subsets 
S C {1, 8}, creates the submatrix that consists of the columns of A given by the 
numbers in S, and computes the determinant of this submatrix. 

2. Construction of the Jacobian ideal 

In order to carry out explicit computations in the cohomology of X one replaces the 
intersection of four quadrics given by the matrix A by a certain toric hypersurface 
with defining equation 

F = yifi + V2f2 + 2/3/3 + V4f4 with f 1 ,...,f 4 e K[x ,...,x 7 }- 

as outlined in section 4 of [T]. The next step is to define the polynomial ring and 
the Jacobian ideal associated to this hypersurface. The polynomial ring is created 

by 

PR : =PolynomialRing (K , 12 , "glex" ) ; 
xn: = [] ; 

for i:=0 to 7 do 

xn[i+l]:="x" cat IntegerToString(i) ; 
end for; 
yn: = []; 

for i:=l to 4 do 

yn[i]:="y" cat IntegerToString(i) ; 
end for; 

AssignNames(~PR,xn cat yn) ; 
x:=function(i) 

return PR. (1+i) ; 
end function; 
y :=function(i) 

return PR. (8+i) ; 
end function; 

The polynomial ring PR has 12 variables xq, ...,x?,yi, ...,2/4 which are accessed by 
x(i) and y(j) and diplayed as xi, yj, respectively. Now the defining equation of 
the toric hypersurface can be computed from the coefficient matrix A. 



f : = □ ; 

for i:=l to 4 do 
f [i] :=PR!0; 
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for j : =0 to 7 do 

f [i] :=f [i]+A[i] [j+l]*x(j)-2; 

end for; 
end for; 
F:=PR!0; 

for i:=l to 4 do 

F:=F+y(i)*f [i] ; 
end for; 

The resulting equation is now stored in the variable F. Finally, we define the Jaco- 
bian ideal associated to this equation. 

I: = []; 

for i:=l to 12 do 

I : =Append(I .Derivative (F, i) ) ; 
end for; 
I:=Ideal(I) ; 
I :=GroebnerBasis(I) ; 
I:=IdealWithFixedBasis(I) ; 



3. Deducing a cohomology basis 

In this section we describe how to determine a basis for the cohomology of X as 
desribed in section 4 of pQ . In section 2 we defined an action on the ring PR by a 
certain finite group Ni . The cohomology of X is given by the iVi-fixed part of this 
ring modulo the Jacobian ideal. It is a graded iC-algebra 

R = @R P 

where the component R p is generated by the monomials with total degree 2p in the 
variables Xo, ...,xj and total degree p in yi, j/4. It follows that a basis of R can 
be computed by the following procedure: 

(i) Enumerate all monomials of x-degree 2p and y-degree p. 

(ii) Discard all monomials which do not agree with their normal form with 
respect to the Jacobian ideal. This yields a basis of R p . 

(iii) Discard all monomials which are not fixed by the action of N±. The re- 
maining monomials constitute a basis of R p . 

For (i) we define two additional polynomial rings with 8 and 4 indeterminates, 
respectively. 

XR:=PolynomialRing(K,8) ; 
YR:=PolynomialRing(K,4) ; 

We also provide functions which map the indeterminates of XR to xq, X7 and 
the indeterminates YR to y\, j/4, respectively. Using these functions, we are now 
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able to enumerate all monomials of a;-degree i and y-degree j as follows: First, we 
use the Magma function MonomialsOf Degree in order to produce all monomials 
of degree i in XR and degree j in YR. The functions below cast these monomials 
into elements of PR, and the desired monomial set is obtained by multiplying any 
element of the first set with any element of the second. 

xmon2mon : =f unction(xmon) 
res:=PR!l; 
for i:=0 to 7 do 

exp : =Degree(xmon, i+1) ; 
res : =res* (PR. (i+1) ) "exp; 
end for; 
return res; 
end function; 

ymon2mon : =f unction(ymon) 
res:=PR!l; 
for i:=l to 4 do 

exp : =Degr e e ( ymon , i ) ; 
res:=res*(PR. (i+8))~exp; 
end for; 
return res; 
end function; 

Step (ii) is based on the following simple function. 

isbasiselt : =f unction (mon) 

nf : =NormalForm (mon , I ) ; 

return (nf eq mon) ; 
end function; 

Finally, step (iii) is carried out by the function below. It makes use of the fact that 
the group Ni is generated by a a , where a runs through the set 

{ei + e 2 , e 2 + e 3 , e 7 + e 8 , e 8 + ei} 

(see section 4 of [T] for the notation). 

isHinvar : =f unction (mon) 
cond:=true; 
for i:=0 to 6 do 

degl :=Degree (mon, i+1) ; 

deg2 : =Degree (mon , i+2) ; 

cond:=cond and IsEven(degl+deg2) ; 
end for; 

degl : =Degree (mon, 1) ; 
deg2 : =Degree (mon , 8) ; 
cond:=cond and IsEven(degl+deg2) ; 
return cond; 
end function; 
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Based on these functions, the following loop now determines a basis for every com- 
ponent the graded ring. The result is the nested array basis in which basis lp+ 1] 
contains the basis elements of R p , given as elements of the ring PR. 

basis : = [] ; 

for i:=0 to 2 do 

basis [i+1] : = [] ; 

xdeg : =2*i ; ydeg : =i ; 

xmons : =MonomialsOf Degree (XR.xdeg) ; 

ymons : =MonomialsOf Degree (YR,ydeg) ; 

for xmon in xmons do 

for ymon in ymons do 
mon : =xmon2mon (xmon) *ymon2mon (ymon) ; 
cond : =isHinvar (mon) ; 
cond:=cond and isbasiselt (mon) ; 
if cond then 

basis [i+1] : =Append(basis [i+1] ,mon) ; 
end if ; 
end for; 
end for; 
end for; 

basis[4] : = [x(7) ~6*y (4) ~3] ; 

In order to save time, we define the component basis [p + 1] directly instead of 
using the procedure. The reason is that the space of polynomials of bidegree (6, 3) 
is already quite large, so the computation would take several seconds. 

4. The first characteristic subvariety 

We briefly recall the definition of the characteristic subvariety IZi from section 3 
in pQ. Let ei, ...,eg denote a basis of R\ and e[, ...,e 9 a basis of i?2- Taking the 
multiplication map to its dual yields a linear map 

H* :R* 2 — ► S 2 (R\). 

Fix a bijection ip : {(i,j) G IN 2 | 1 < i < j < 9} — ► {1, ...,45}. Furthermore, let 
A 6 x 9x45 denote a representation matrix of /i* with respect to (e^)*, (eg)* and 
b(p(i,j) '■= e i e j- Then IZi is the subvariety in P 8 given by the equations 

9 9 

ft ■= ^2^2 a £M*j) Z * Z 3 G K l Z l' ■■■' z 9]- 
i=l j=i 

Thus in order to compute the variety, we have to carry out the following steps: 

(i) provide functions for the bijection ip and its inverse ip~ l 

(ii) compute a representation matrix of the multiplication map /i : S 2 (R\) — > 
i?2 with respect to the basis b v u,j) ■— and e[, e 9 

(iii) obtain a representation matrix with respect to (e^)*, (e 9 )* and b* ^ ^ = 
(eiej)* by transposition 

(iv) deduce matrix A by base change from (eiej)* to e*e* 
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(v) write down the defining equations of IZi 
The bijections <p and ip -1 are provided by the following two functions. 

r :=PolynomialRing(GF(2) ,9) ; 
symm2order : =MonomialsOf Degree (r , 2) ; 
symm2order : =SetToSequence (symm2order) ; 
Symm2Spc : =VectorSpace (K,#symm2order) ; 

ijpos:=function(i, j) 
if i gt j then 
i_:=j; j_:=i; 

else 

i_:=i; j_:=j; 
end if ; 

p:=9 - i_ + 1; 
offs:=45 - l/2*p*(p+l) ; 
offs:=offs + (j_-i_) + 1; 
offs:=Integers() loffs; 
return offs; 
end function; 

pos2ij :=function(pos) 

mon: =symm2order [pos] ; 
res : = [] ; 

for i:=l to 9 do 

if Degree (mon, i) eq 2 then 

return [i,i] ; 
end if ; 

if Degree (mon, i) eq 1 then 

res : =Append(res , i) ; 
end if ; 
end for; 
return res; 
end function; 

In order to carry out (ii) we need a function which computes, for an arbitrary 
polynomial g in the polynomial ring PR, a coordinate vector of length 1 + 9 + 9 + 1 = 
20 of the cohomology class represented by g with respect to basis. 

CSpc : =VectorSpace (K , 20) ; 
poly2vec : =f unction(poly) 

res:=CSpc!0; 

nf : =NormalForm (poly , I ) ; 

offs:=0; 

for p:=0 to 3 do 

for i:=l to #basis[p+l] do 

coe:=MonomialCoefficient(nf .basis [p+1] [i] ) ; 

res[offs+i] :=coe; 
end for; 

of f s :=of f s+#basis [p+1] ; 
end for; 
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return res; 
end function; 

The representation matrix is now obtained by computing the products of all pairs 
of elements in basis [2] . The 45 rows of the matrix are given by the i? 2 -parts of 
the coordinate vectors of these products. 

M: = [] ; 

for 1:=1 to 45 do 
M[l]: = []; 

i j : =pos2i j (1) ; i : =i j [1] ; j : =i j [2] ; 

monl :=basis [2] [i] ; mon2 : =basis [2] [j] ; 

vec : =poly2vec (monl*mon2) ; 

for k:=l to 9 do 

M[l] [k] :=vec[10+k] ; 

end for; 
end for; 
M:=Matrix(K,M) ; 

The implementation of the steps (iii) and (iv) is straightforward, provided one takes 
into account that (e^)* = e*e* and (e^)* = 2e*e*. The final result is stored in 
the matrix named C. 

N: = [] ; 

for 1:=1 to 45 do 
N[l] : = [] ; 
ij :=pos2ij (1) ; 
S:=SequenceToSet(ij) ; 
if #S eq 2 then 
la:=2; 

else 

la:=l; 
end if ; 

for k:=l to 9 do 

N [1] [k] : =la*M [1] [k] ; 

end for; 
end for; 
N:=Matrix(K,N) ; 
C:=Transpose(N) ; 

Finally, the matrix C is used to carry out (v). The resulting projective scheme is 
stored in the variable charvar. 

CSR:=PolynomialRing(K,9) ; 
chareqs : = [] ; 
for 1:=1 to 9 do 

f :=CSR!0; 

for k:=l to 45 do 

i j : =pos2i j (k) ; i : =i j [1] ; j : =i j [2] ; 
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mon:=CSR.i*CSR. j ; 

f :=f+C[l] [k] *mon; 
end for; 
chareqs [1] : =f ; 
end for; 

P8 : =Pro j ectiveSpace (CSR) ; 
charvar : =Scheme (P8 , chareqs) ; 

The dimension and arithemtic genus are obtained by the function calls 
Dimension(charvar) and ArithmeticGenus (charvar) . 

5. The second characteristic subvariety 

The computation of the characteristic subvariety IZ2 works essentially in the same 
way as in the previous section. As before, let e\, ...,eg denote a basis of R\. By 
e[ we denote a vector which spans the 1-dimensional space R3. The dual of the 
multiplication map yields 

\i* : i?g — > S Z (R\). 

In order to write down a representation matrix, we fix a bijection 

: k) e IN 3 I l<i<j'<fc<9} — > {1, 210}. 

Let A £ x lx165 a representation matrix of fi* with respect to e[ and b ip ( i j_k) '■— 
e*e*e* k . Then H2 is by definition a hypersurface in P 8 given by the quation 

999 

• 9:= EEE a l, V {i,j,k)ZiZjZ k e K\z x , z 9 }. 
i=l j=i k=j 

Thus we can proceed as follows. 

(i) provide functions for the bijection tp and its inverse (p^ 1 

(ii) compute a representation matrix of the multiplication map \x : S 3 (Ri) — > 
i?3 with respect to the basis & ¥ >(i,j,fe) : = ZiZjZk and e[ 

(iii) obtain a representation matrix with respect to (e^)* and 6* ^ . ^ — (eiejek)* 
by transposition 

(iv) deduce matrix A by base change from (eiejek)* to e*e*e*. 

(v) write down the defining equation of IZ2 

Step (i) is achieved by 

r:=PolynomialRing(GF(2) ,9) ; 
symm3order : =MonomialsOf Degree (r , 3) ; 
symm3order : =SetToSequence (symm3order) ; 
Symm3Spc : =VectorSpace (K,#symm3order) ; 

ijkpos : =f unction (i , j ,k) 
S:=Sort([i, j ,k] ) ; 
i_:=S[l] ;j_:=S[2] ;k_:=S[3] ; 
offs:=0; 

for 1:=1 to i_-l do 
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p:=10 - 1; 

offs:=offs + l/2*p*(p+l); 

end for; 

p0:=9-i_+l; 

j0:=j_-i_+l; 

k0:=k_-i_+l; 

offs2:=0; 

p0_:=p0 - jO + 1; 

offs2:=l/2*p0*(p0+l) - l/2*p0_*(p0_+l) ; 

offs2:=offs2 + (kO - jO) + 1; 

off s:=off s+offs2; 

off s:=Integers() !offs; 

return offs; 

end function; 

pos2ijk:=function(pos) 
mon : =symm3order [pos] ; 
res: = [] ; 

for i:=l to 9 do 

if Degree (mon, i) eq 3 then 

return [i , i , i] ; 

end if ; 

if Degree (mon, i) eq 2 then 
res:=res cat [i,i]; 
end if ; 

if Degree (mon, i) eq 1 then 

res:=Append(res,i) ; 

end if ; 

end for; 

return res; 

end function; 

For step (ii) we can use the function poly2vec from the previous section. The 
representation matrix of the multiplication map \x is then computed by 

M: = [] ; 

for 1:=1 to 165 do 
M[l] : = [] ; 

ijk:=pos2ijk(l) ;i:=ijk[l] ; j :=ijk[2] ;k:=ijk[3] ; 

monl : =basis [2] [i] ; mon2 :=basis [2] [j] ; 

mon3 : =basis [2] [k] ; 

vec :=poly2vec (monl*mon2*mon3) ; 

M[l] [1] :=vec[20] ; 

end for; 

M:=Matrix(K,M) ; 

For step (iii) and (iv), notice that (e 3 )* = (e*) 3 , (efej)* = 3(e*) 2 e* and (eiejek)* = 
6e*e*e£ for k pairwise distinct. 

N: = []; 

for 1:=1 to 165 do 
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N[l] : = [] ; 

ijk:=pos2ijk(l) ; 

S:=SequenceToSet(ijk) ; 

case #S: 

when 1 : 

la:=l; 

when 2 : 

la:=3; 

when 3 : 

la:=6; 

end case; 

N[l] [1] :=la*M[l] [1] ; 
end for; 
N:=Matrix(K,N) ; 
C:=Transpose(N) ; 

Now the characteristic variety is obtained by 

CSR:=PolynomialRing(K,9) ; 

chareq:=CSR!0; 

for 1:=1 to 165 do 

i jk : =pos2ijk(l) ; 

i:=ijk[l] ;j:=ijk[2] ;k:=ijk[3] ; 

mon : =CSR . i*CSR . j *CSR . k ; 

chareq:=chareq+C [1] [l]*mon; 
end for; 

P8 : =Pro j ectiveSpace (CSR) ; 
charvar : =Scheme (P8 , chareq) ; 



6. Computation of the Higgs field 

For the implementation of the Plethysm method (see section 3 in [I]) it is neccessary 
to compute a representation matrices A w £ x 20x20 of the maps 

fi w : R — ► R, x i— » wx for w G R\ 

with respect to the basis computed in section[3l The idea is simple: Multiply every 
basis element with some fixed w € R±, use the function poly2vec from section [4] 
in order to transform the products into row vectors and form a matrix out of these 
rows. This task is carried out by the following code. The result is an array called 
thetamats whose elements are the representation matrices corresponding to the 9 
basis vectors of R\ . 

thetamats : = [] ; 
for j : =1 to 9 do 

the:=basis [2] [j] ; 
thetamat : = [] ; 
for p:=0 to 3 do 

for i:=l to #basis[p+l] do 
poly : =the*basis [p+1] [i] ; 
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vec : =poly2vec (poly) ; 
vec :=ElementToSequence(vec) ; 
thetamat :=Append(thetamat, vec) ; 
end for; 
end for; 

thetamat :=Matrix(K,thetamat) ; 
thetamats : =Append(thetamats , thetamat) ; 
delete thetamat; 
end for; 



7. The induced Higgs field on symmetric 2-space 

Each map fi w : R — > R induces an endomorphism on S 2 (R). The aim of this section 
is to represent the elements of S 2 (R) by vectors in K 210 , and to compute the action 
of fi w with respect to this representation. For the first part, we only have to fix a 
bijection 

^{^OelN 2 | l<i<j<20}— ►{1,2,...,210}. 
If 61, ...,620 denotes a basis of R, then the isomorphism S 2 (R) — > K 210 is given by 
bibj 1 ► e<p(%,j)- The following two functions provide such a bijection. 

r :=PolynomialRing(GF(2) ,20) ; 
symm2order : =MonomialsOf Degree (r , 2) ; 
symm2order : =SetToSequence (symm2order) ; 
Symm2Spc : =VectorSpace (K,#symm2order) ; 

ijpos : =function(i , j ) 
if i gt j then 
i_:=j; j_:=i; 

else 

i_:=i; j_:=j; 
end if ; 

p:=20 - i_ + 1; 
offs:=210 - l/2*p*(p+l); 
offs:=offs + (j_-i_) + 1; 
offs:=Integers() !offs; 
return offs; 
end function; 

pos2ij : =f unction(pos) 

mon: =symm2order [pos] ; 
res: = [] ; 

for i:=l to 20 do 

if Degree (mon, i) eq 2 then 

return [i , i] ; 
end if ; 

if Degree (mon, i) eq 1 then 
res:=Append(res,i) ; 

end if ; 
end for; 
return res; 
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end function; 

Given an basis vector w € R\ and some 6 K 210 = S 2 (R), our next task is 
two compute the image n w (ek)- By the above isomorphism, corresponds to an 
element of the form bibj, and the induced map is given by 

H{bibj) = bifi(bj) + fjb(bi)bj. 

It means that we have to proceed as follows: 

(i) Use the matrices in thetamats in order to compute the images n(bi) and 
^(bj) as elements in if 20 . 

(ii) Compute the product bi[i(bj); if [i(bj) = Y^k=i a kbk, then 

20 



fe=i 

Use the bijection ip in order to map this element to K 210 . 
(iii) Do the same for /i(bi)bj, and return the sum. 
This is carried out by the following Magma code. 

symm2imthe : =f unction(th,pos) 

i j : =pos2i j (pos) ; i : =i j [1] ; j : =i j [2] ; 

imi : =CSpc ! ; imi [i] : =1 ; 

imj : =CSpc ! ; imj[j]:=l; 

imj :=imj*thetamats [th] ; 

rest :=Symm2Spc ! 0; 

for k:=l to 20 do 

if not (imj [k] eq 0) then 

rest [ijpos(i,k)] :=imi[i]*imj [k] ; 

end if ; 
end for; 

imi : =CSpc ! ; imi[i]:=l; 
imi : =imi*thetamats [th] ; 
imj : =CSpc ! ; imj[j]:=l; 
res2:=Symm2Spc!0; 
for k:=l to 20 do 

if not (imi [k] eq 0) then 



end if ; 
end for; 
res:=resl+res2; 
return res; 
end function; 

By linearity we can now compute the image of arbitrary elements in K 210 . 

symm2imthe_ :=function(th,vec) 
res : =Symm2Spc ! ; 
seq: =ElementToSequence (vec) ; 
for i:=l to #seq do 




res2 [i jpos (k, j )] :=imi[k]*imj [j] ; 
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imvec : =symm2imthe(th, i) ; 

res : =res+vec [i] *imvec ; 
end for; 
return res; 
end function; 

For the application of the plethysm method, it is essential to carry out the following 
task: Given a subspace U < S 2 (R), compute the sum 

(U) © ii W2 (U) © • • • © fj, W9 (U) 

of the images of all Higgs field maps. (Here w\, wg denotes a basis of the graded 
subspace R\.) We use the built-in linear algebra functions of Magma in order to 
provide this. 

symm2image : =f unction (U) 
B:=Basis(U) ; 
ims : = [] ; 

for bvec in B do 

for th:=l to 9 do 

ims : =Append ( ims , symm2imthe_ (th , bvec) ) ; 

end for; 
end for; 

return sub<Symm2Spc I ims> ; 
end function; 

8. The graded subspaces of symmetric 2-space 
The grading R — (B p=0 R P on the Jacobian ring induces a natural grading 

S 2 (R) = Y,S 2 (R) P 

on symmetric 2-space. For the plethysm method, we have to compute the iterated 
images of the graded subspaces S 2 (R) P under the Higgs field elements. In particular, 
for < p < 6 we have to determine the image of S 2 (R) P under the isomorphism 
S 2 (R) P = K 2W . We provide a function which returns for each p a set {ii, ...,ij p } 
such that the image of S 2 (R) P in K 2W is spanned by the basis vectors e^, ...,e^ . 

It seems to be the easiest strategy to compute the index set for each p "by hand" . 
Let b\ 1 620 denote a basis of R. We assume that the indices are chosen such that 
Ro is spanned by bi, R\ is spanned by b2,---,bio, i?2 by the elements bu,...,b\g 
and finally R 3 by 6 2 o- We first consider p = 0. In this case S 2 (R)o is generated by 
the image of i?o © Ro, which means that the image of S 2 (R)o in K 210 is spanned 
by e v (i .!). The space S 2 (R)i is generated by Ro © so the image is spanned 
by b\bi — e v n t i) with 2 < i < 9. For p = 2 we have to consider the image of 
(Ro © R2) © (Ri © Ri) etc. If we work through all degrees up to 6, we end up with 
the following function. 



DISPROOF OF MODULARITY 



symm2degind : =f unction (p) 
case p: 

when : 

return [ijpos (1 , 1)] ; 
when 1 : 

res : = [] ; 

for i:=l to 9 do 

res:=Append(res,ijpos(l,l+i)) ; 
end for; 
return res ; 
when 2 : 

res: = [] ; 

for i:=l to 9 do 

res:=Append(res,ijpos(l,10+i)) ; 
end for; 

for i:=l to 9 do 

for j :=i to 9 do 

res : =Append (res , ijpos ( 1+i , 1+ j ) ) ; 
end for; 
end for; 
return res ; 
when 3: 

res:=[ijpos(l,20)] ; 
for i:=l to 9 do 

for j :=1 to 9 do 

res : =Append (res , ijpos ( 1+i , 10+ j ) ) ; 
end for; 
end for; 
return res ; 
when 4: 

res: = [] ; 

for i:=l to 9 do 

res : =Append(res , ijpos (1+i ,20) ) ; 
end for; 

for i:=l to 9 do 

for j :=i to 9 do 

res : =Append (res , ijpos ( 10+i , 10+ j ) ) ; 
end for; 
end for; 
return res ; 
when 5 : 

res : = [] ; 

for i:=l to 9 do 

res : =Append(res , ijpos (10+i , 20) ) ; 
end for; 
return res ; 
when 6 : 

res:=[ijpos(20,20)] ; 
return res; 
end case ; 
return -1; 
end function; 
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The indices provided by symm2degind now admit the computation of the p-th 
graded subspace of S 2 (R). 

symm2degspc : =function(p) 
basis : = [] ; 

inds : =symm2degind (p) ; 
for i:=l to #inds do 

vec : =Symm2Spc ! ; 

vec [inds [i] ] : =1 ; 

basis :=Append(basis, vec) ; 
end for; 

return sub<Symm2Spc I basis> ; 
end function; 

The above function now have to be invoked in the following way in order to make 
the plethysm method work. 

U5 1 : =symm2degspc ( 1 ) ; 
U42:=symm2image(U51) ; 
U33:=symm2image(U42) ; 
Dimension (U33) ; 

The result is 78, whereas the modularity assumption would have implied a dimen- 
sion < 65. 
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